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Abstract 

Topological defects in ordered states with spontaneously broken symmetry often have unusual 
physical properties, such as fractional electric charge or a quantised magnetic field-flux, originating 
from their non-trivial topology. Coupled topological defects in systems with several coexisting or- 
ders give rise to unconventional functionalities, such as the electric-field control of magnetisation in 
multiferroics resulting from the coupling between the ferroelectric and ferromagnetic domain walls. 
Hexagonal manganites provide an extra degree of freedom: In these materials, both ferroelectricity 
and magnetism are coupled to an additional, non-ferroelectric structural order parameter. Here 
we present a theoretical study of topological defects in hexagonal manganites based on Landau 
theory with parameters determined from first-principles calculations. We explain the observed flip 
of electric polarisation at the boundaries of structural domains, the origin of the observed discrete 
vortices, and the clamping between ferroelectric and antiferromagnetic domain walls. We show that 
structural vortices induce magnetic ones and that, consistent with a recent experimental report, 
ferroelectric domain walls can carry a magnetic moment. 



Recent experimental and theoretical advances in the field of multiferroics have clarified 
the microscopic mechanisms of coupling between ferroelectricity and magnetism in bulk ma- 
terials and led to the discovery of many families of materials in which ferroelectricity is 
induced by a spin ordering. Examples include the orthorhombic rare-earth manganites, 
spinels, hexaferrites and delafossites, which usually have spin orderings of the cycloidal or 
conical spiral typeP^ The resulting electric polarisation is highly susceptible to an applied 
magnetic field and can be easily rotated or reversed However, this magnetically-induced 
polarisation is usually too small to allow manipulation of spin states by an applied voltage. 
Much larger electric polarisations are found in multiferroics such as BiFeC>3 and the hexago- 
nal rare-earth manganites, in which ferroelectricity results not from spin ordering but from 
electronic and lattice instabilities.^! Yet in these materials, the electric control of magnetism 
is not straightforward,^^] s i nce the direction of spins in the magnetically ordered state is 
not correlated with the sign of the macroscopic electric polarisation.^ 

While enhancing bulk couplings between polarisation and magnetism is difficult, practical 
switching of a ferroic order parameter with an applied field invariably involves motion of 
the domain walls. Magnetoelectric switching therefore depends crucially on interactions 
between ferroelectric and ferromagnetic domain walls, which are not as well understood. 
In this context, the observed clamping between ferroelectric and antiferromagnetic domain 
walls in multiferroic hexagonal manganites 14 provides a unique prototype for investigation. 

The hexagonal manganites, i?Mn03, where R denotes a small-radius rare earth ion (Dy, 
Ho, Er, Tm, Yb, Lu), Y or Sc, are improper ferroelectrics; electric polarisation appears as 
a by-product of a primary structural transition.^ The crystal structure consists of corner- 
sharing MnOs trigonal bipyramids, which form triangular layers, separated by layers of R 
ions. The structural transition above 1000 K results in periodic tilts of the MnOs bipyra- 
mids and displacements of the R ions along the c axis normal to the layers.^^ This periodic 
v3 x lattice distortion makes the size of the unit cell three times larger and is often 
referred to as the trimerisation transition (see Fig. [IJa-c)). The anharmonic coupling be- 
tween the trimerisation mode and a polar optical phonon mode induces the observed electric 
polarisation along the c axis, P c ~ 5.5 /xC ■ cm" 2 P 

At much lower temperatures, ~ 100 K, an antiferromagnetic ordering of Mn spins emerges. 
While there is a large body of evidence for the strong interplay between the microscopic 
spin, charge and lattice degrees of freedom in hexagonal manganites] 16 * 18 ^ the sign of the 
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overall antiferromagnetic order parameter is decoupled from the direction of the macroscopic 
electric polarisation, as such correlation is forbidden by symmetry. It therefore came as a 
surprise when non-linear optical measurements of YMnC>3 demonstrated that ferroelectric 
domain walls are locked to magnetic ones.^ Furthermore, this clamping was found to be 
non-reciprocal, as "free" magnetic domain walls, not associated with the electric polarisation 
reversals, were also observed. Proposals for the clamping mechanism have included strain 
mediation^ and renormalisation of spin interactions at ferroelectric walls .122123] 

Recently, new observational evidence has shed light on the mechanism for clamping. A 
combination of conducting atomic-force microscopy and transmission-electron microscopy 
demonstrated that the ferroelectric domain walls are pinned to the boundaries of the struc- 
tural domains that appear upon transition to the trimerised statel^^ES These measure- 
ments also revealed intricate patterns of unusual line defects. These so-called 'cloverleaf 
defects,^^ at which six different structural and ferroelectric domains merge, have also been 
seen in piezoresponse force microscopyP^H In Reff2S1 it was suggested that the line defects 
are discrete analogues of vortices, and that the change of polarisation sign at structural 
domain boundaries is a consequence of the special form of the coupling between the lattice 
distortion and electric polarisation originating from the "geometric" nature of ferroelectric- 
ity in hexagonal manganites.^^ These results put the clamping between the ferroelectric 
and antiferromagnetic domain walls into an entirely new perspective. 

In this paper we study the interplay between structural, ferroelectric and magnetic defects 
in hexagonal manganites using an expansion of the free energy in powers of the correspond- 
ing order parameters and their gradients. Based on symmetries of ordered states of these 
materials, this approach allows us to identify stable topological defects and describe their 
mutual interactions in the most economical way. In ReflTTl Fennie and Rabe discussed the 
Landau theory of improper ferroelectricity for spatially uniform states of hexagonal man- 
ganites. Extracting parameters of the Landau expansion from first-principles studies of 
YMnC>3, they showed that polarisation emerges due to a non-linear coupling to the trimeri- 
sation mode. We extend this theory to inhomogeneous topological defects and include spin 
degrees of freedom to study effects of the structural domain walls and vortices on magnetic 
ordering. 

In the trimerised state three neighbouring MnOs bipyramids tilt towards (or away from) 
their common equatorial oxygen atom (see Fig. ma-c))PAs a consequence of the hexago- 
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FIG. 1: Structural and magnetic angles. a,b,c, Projections of the displacements of apical 
oxygen ions on the ab plane in the trimerised state indicated by arrows. One ion (blue arrow) is 
chosen to define the trimerisation phase <£. Shown are the a+ state with $ = (a), 7_ state with 
$ = 7r/3 (b) and the (3+ state with <3? = ^ (c). Also shown are the displacements of apical oxygen 
ions in an adjacent Mn-0 layer (lighter triangles). d,e,f,g, The four magnetic states of hexagonal 
manganites with the spin directions indicated by red arrows and the corresponding values of the 
angles (xi,X2), where xi,2 = ipi,2 — $ (see text), h, The angles tpi and V>2 describing the rotations 
of spins in magnetic domain walls (ipi,ip2 > correspond to the clockwise(anticlockwise) rotation 
in even(odd) Mn layers). 

nal structure of Mn-0 layers, there are six distinct trimerised states, corresponding to six 
degenerate minima of the lattice energy. Being a periodic lattice modulation in a layered 
system, the trimerisation is described entirely by the amplitude Q and phase The phys- 
ical meaning of the phase $ is the azimuthal angle describing the in-plane displacements 
of apical oxygens (see Fig. [I|a-c)). The minimal-energy states can then be labelled by the 
six values of the phase: 0, ±|,±^ and ir. At the structural domain boundaries $ varies 
spatially between two of these six values. 

Microscopically, the trimerisation is the condensation of the zone-boundary K% mode with 
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wave vector q = (1/3, 1/3, 0), which breaks the P6s/mmc symmetry of the undistorted phase 
lowering it to YQ^cm. Similarly, the spontaneous electric polarisation P c is proportional to 
the amplitude of the zone-centre mode V with symmetry. This polar mode is stable in 
the PQ^/mmc structure, but is non-linearly coupled to the unstable K% mode, and therefore 
appears together with the trimerisation. 
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TABLE I: Transformations of the trimerisation phase the spin angles tpi and "02 > the electric 
polarization P c and the magnetic field H c under the generators of the VQ^/mmc space group 
describing the high-temperature phase: translation S a = (x + 1, y, z), three-fold axis 3 C = (— y, x — 
y,z), two-fold screw axis 2 C = (— x, — y, z + 1/2), mirror plane m a+ b = (—y,—x,z), inversion 
I = (—x, —y, —z), and the time reversal operation T. 

The free-energy expansion in powers of Q, V, and their gradients, 

/= ^ 2 + ^ 4 + ^(c + c'cos6$) 

- g Q*V cos 3$ + 9 -Q 2 V 2 + ^V 2 (1) 
+ \ i s Q ( d iQ d iQ + Cfd&d&) + SrdpPdpP] , 

i=x,y,z 

is obtained using the transformation properties of the trimerisation phase $ and the polari- 
sation P c under the generators of the high-temperature space group summarised in Table |T} 
We consider only the lowest-order stiffness terms accounting for the energy cost of spatial 
variations of Q and V. (x, y) are the Cartesian coordinates in the ab plane (see Fig. [T]), and 
by symmetry, Sq = Sq and Sp = s y v . 

The trimerisation phase $ and the stiffness terms, sq and s-p, not considered by Fennie 
and Rabep^ play an important role in the theory of topological defects. In particular, 
the form of the non- linear coupling, — gQ 3 V cos 3$, giving rise to improper ferroelectricity, 
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implies that for g > the electric polarisation induced in the states with $ = 0, +-^- and 
— y is positive (the a> + , /3 + and 7+ phases) while for +f,7T and — | it is negative (the 
7_, a_ and /3_ phases). In other words, neighbouring trimerisation phases, separated by 
A$ = |, have opposite electric polarisations. 

a = -2.626 eV- A -2 b = 3.375 eV • A" 4 c = 0.117 eV • A~ 6 d = 0.108 eV • A - ' 
a P = 0.866 eV-A -2 p = 1.945 eV • A" 4 g' = 9.931 eV • A -4 
s z Q = 15.40 eV s x Q = 5.14 eV sp = 52.70 eV sp = -8.88 eV 

TABLE II: Parameters of the phenomenological expansion of the free energy Eq. ([2]) obtained from 
ab initio calculations. All parameters are calculated per unit cell of the trimerised lattice with the 
volume V = 365.14 A 3 containing 6 formula units. The polarization P c is related to the amplitude 
of the polar mode V by P c = V~ l Z*V, where Z* = 9.031e is the effective charge of the polar 



mode. Because of the negative stiffness, sp, of the polar mode, the term \tp 
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with tp = 73.56eV • A was added to Eq.(l) in order to calculate the structure of the domain wall 
shown in Fig. |3^a). 

The numerical values of the parameters a, b, etc. in Eq. ^ for YMnC>3 are listed in 
Table [TT} In order to determine them, we performed ab initio supercell calculations for the 
various lattice distortions (see Appendix A). In particular, Fig. [2]^a) shows the variation of 
the total energy with respect to the magnitude of Q for two chosen values of $ and V = 0. 
It is clear that in the physically relevant range of Q, the total energy depends very weakly 
on $. The ^-dependence of the energy f(Q, $, V — 0) first appears in the sixth-order of the 
expansion in powers of Q and the corresponding coefficient d is small. The energy landscape 
in the (Q, $) plane therefore essentially has a rotationally invariant Mexican Hat shape with 
no barriers separating the six structural domains. 

In reality, the ^-rotation of the tilted bipyramids is not a zero mode, because an addi- 
tional cos 6$ term is generated by minimising / with respect to V and eliminating V from 
Eq.(2), which lowers the energy by 9 ||^p^y^ • Fig- 2^ c ) shows the dependence of the total 
energy on V for $ = (at which the energy has minimum) and $ = 7r/6 (at the top of the 
barrier separating two minimal-energy states). Clearly, the anharmonic coupling between 
Q and V results in a strong $ dependence of the energy. Therefore, even though the emer- 
gence of ferroelectricity in hexagonal manganites is improper, the coupling of trimerisation 
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FIG. 2: Extraction of model parameters from ab initio calculations, a, Variation of 
total energy with respect to K3 mode amplitude for two trimerisation angles, = and <I> = 
Parameters a, b, c and c' are extracted, and owing to the weak angular dependence of the energy d is 
small, b, Variation of total energy about the high-symmetry P63/mmc structure with respect to the 
amplitude of the (polar) mode, V. The polar mode is stable and does not spontaneously emerge 
at T = OK. c, Coupling between K% and T2- modes for two different trimerisation angles. For 
trimerisation angles of 0, ±|, ±^ and ir the anharmonic coupling leads to a non-zero polarisation 
and a total energy lowering of ~ 26meV per formula unit. On the other hand, intermediate 
trimerisation angles do not allow polarisation to develop. The g and g 1 parameters are extracted 
from these data, d, Contour plot of the free energy of uniformly trimerised states as a function of 
Q and Here, P z has been optimized for each Q, The trajectory Q(&) (white dashed line) 
connecting two neighbouring energy minima corresponds to the lowest-energy structural domain 
wall. 

to polarisation is the only factor that determines the energetic barriers between different 
trimerised states, and is responsible for replacing the accidental continuous XY symmetry 
of / with the discrete Z$ symmetry. 

This has a strong effect on the structure of topological defects in the trimerised state, 
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FIG. 3: Structural topological defects. a,b, The variation of the trimerisation angle $ (mea- 
sured in radians, thick blue line), trimerisation amplitude Q (measured in A, thin green line) and 
the polar mode V (measured in A, dashed red line) across the lowest-energy domain wall normal to 
the ab plane (a) and parallel to the ab plane (b). c, The vortex-antivortex pair. In the structural 
vortex/antivortex the angle increases/decreases by 2tt along the loop encircling the vortex core 
in the counter-clockwise direction (the 6 trimerisation states are indicated by colour), d, The 
'topological stripe domain state' in a thin film of ferroelectric hexagonal manganite parallel to the 
ab plane with the alternating polarisation (white symbols) along the c axis. A strain indicated by 
green arrows results in the monotonic increase of the trimerisation angle in the direction normal 
to the stripes. 
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which can be described as trajectories in the (Q, $)-plane minimising the energy for given 
initial and final conditions. For example, a structural domain wall corresponds to a path 
connecting two energy minima. The shortest path connecting two 'neighbouring' states 
whose trimerisation angles differ by A$ = ±7r/3 (see Fig. [2|d)) is the lowest-energy domain 
wall. This path follows the bottom fold of the Mexican Hat where the potential barrier 
between the two minima is the lowest, so that the amplitude of the trimerisation Q in the 
wall is close to its bulk value. Figures |3^a-b) show the coordinate dependence of $ and 
Q across domain walls obtained by numerical free-energy minimisation. The domain wall 
width is 5 — 10 A and Q at the domain wall is reduced by about 10%. 

Since the neighbouring energy minima separated by A$ = ±7r/3 have anti-parallel elec- 
tric polarisations, the structural domain wall is at the same time a ferroelectric domain wall 
(see Fig. [3]^a-b)). The improper nature of ferroelectricity in hexagonal manganites forbids 
purely ferroelectric domain walls, i.e. the reversals of P c within one structural domain, since 
the sign of P c is uniquely determined by the sign of cos 3$. Furthermore, the structural 
domain walls with A$ = 27r/3, separating states with the same electric polarisation, are 
unstable: they 'decay' into lowest-energy domain walls with A$ = 7r/3 which we find to 
repel each other when they overlap. This explains the clamping between the structural and 
ferroelectric domain walls observed by Choi et al^ 

Spherical or cylindrical structural domains are unstable against shrinking, which explains 
the scarcity of closed domain wall lines in hexagonal mang anitesP^ 1 To be stable, the 
domain walls must terminate either at the surface of the sample or at another type of 
stable topological defect - the structural vortex or anti- vortex, shown in Fig. |3](c). At the 
vortex line, where the trimerisation amplitude Q vanishes, all six structural domains meet 
in such an order that the trimerisation phase $ changes by 2n around a contour encircling 
the vortex lineP^ Such a defect cannot be unwound and can only be annihilated by an 
anti- vortex, around which the phase changes by — 2tt (see Fig. [3|c)). Away from the core, 
the trimerisation phase $ varies strongly only at the six radial domain walls. The electric 
polarisation changes sign at each domain wall and varies six times along a loop encircling the 
vortex line. These vortices and anti-vortices are the 'cloverleaf defects' observed in RefJ2H 

Fig. [3|c) shows a vortex-anti-vortex pair configuration obtained by minimising the en- 
ergy for a given distance between these defects. The domain walls diverge radially from the 
vortex/anti- vortex core with the 60°-angle between neighbouring domain walls. Far from 



9 



the core they bend and become parallel to minimise the total length of the structural bound- 
aries, which gives rise to a linear potential between the discrete vortices as opposed to the 
logarithmically growing potential for continuous vortices.^ Despite this confining potential, 
vortex lines and domain walls form dense networks^ 2 -^ that are snapshots of states close to 
critical temperature capturing the formation of these topological defects by large thermal 
fluctuations.™ 11 

A different type of topologically stable domain pattern is uncovered by considering the 
lowest-order coupling of the inhomogeneous trimerisation to strains, which has the form 

/strain = -GQ 2 [(u xx - u yy ) d x <& - 2u xy d y $] , (2) 

where (x, y) are the Cartesian coordinates in the ab plane. With such a coupling, a par- 
allel array of structural domain walls, each with the same increment of the trimerisation 
angle (A<£> = ±7r/3), is topologically stable (see Fig. |3[d)). Because of the alternating 
electric polarisation at the structural walls, such a "$-staircase" is at the same time a fer- 
roelectric stripe domain state, as stabilised in thin films by the long range dipole-dipole 
interactions J 32 * 33 ! These interactions are, however, insensitive to the sign of A<£> at the do- 
main walls, whereas the applied strain selects the direction normal to the walls, in which $ 
increases monotonically. 

Next we address the coupling between the magnetism and the structural/ferroelectric 
domain walls. The spins on the Mn ions in hexagonal manganites order in one of four 
different magnetic states: Ai, A2, Bi, and B2, shown in Fig. []Jd-g).— Their origin can be 
understood by considering the hierarchy of interactions between the magnetic moments on 
Mn sites. By far the strongest is the antiferromagnetic exchange between neighbouring spins 
in the triangular layers of Mn ions, which leads to the non-collinear 120° spin ordering. The 
magnetic easy-plane anisotropy, and the anti-symmetric Dzyaloshinkii-Moriya interaction 
with Dzyaloshinskii vector along the c axis, confine spins to the ab plane. The local in-plane 
anisotropy axes, also favouring the 120° spin angle on neighbouring Mn sites, selects either 
the "radial" (as in the A 2 and B x phases) or the "tangential" (as in Ai and B 2 phases) 
orientation of spins (Fig. [I]). Finally, the interlayer exchange interactions, which are more 
than two orders of magnitude weaker than the intra-layer interactions, lead to either even 
("A" phases) or odd ("B" phases) symmetry under the two-fold screw rotation 2 C P^ 

Due to the strong in-plane exchange, the angle between neighbouring spins remains close 
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FIG. 4: Antiferromagnetic domain walls in hexagonal manganites. a,b, "Free" 180° walls 
within one structural domain. c,d, Domain walls clamped to structural boundaries. Shown is the 
variation of the spin rotation angles tpi and ip2 (solid red and blue lines, respectively), magnetisation 
M c (solid green line, arbitrary units) and the trimerisation phase (dashed blue line) across the 
domain wall. At the structural boundary and the directions of magnetic easy axes on Mn sites 
change by 60°. The directions of spins adjust over much longer distances by rotating over the 60° 
(c) or -120° (d) angle. 

to 120° throughout a magnetic domain wall. Furthermore, to match the directions of spins 
on both sides of the wall with the local in-plane anisotropy axes, the spins in the domain 
wall must rotate in the ab plane. Therefore, the structure of magnetic domain walls can be 
described by two angles, (ipi,ip2), quantifying the rotation of spins around the c axis in the 
even and odd layers respectively, 1 21 1 34 1 as shown in Fig. [TJli) . The transformation properties 
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of ipi and ip2 (see Table |TJ) determine the form of the magnetic free-energy density: 

W^i, fa, H c ) = S [{d^f + («9^ 2 ) 2 ] + A [sin 2 xi + sin 2 X2 ] 

-C+ cos(xi + X2) - C_ cos(xi - X2) - i^M A2 H c (sinxi + sin X2 ) , (3) 

where X i,2 = ^1,2 — *£• The first term originates from the nearest-neighbour exchange. 
Comparing the energies of the 4 uniform phases (see Fig. [l]): fs 2 = —C+—C-, /a 1 = C + +C_, 
/a 2 = 2 A — C + + C_ and = 2A + C + — C_ , we conclude that the A coefficient results from 
the in-plane magnetic anisotropy, C + describes the interlayer exchange interactions, while 
C_ is related to a symmetric exchange anisotropy. The last term in Eq.([3]) describes the 
weak ferromagnetic moment along the c axis in the A 2 phase with {X11X2) = = t( 7r /2, tt/2). 
Terms proportional to spin operators of power higher than 2 are neglected. Equation ^ is 
used to calculate magnetic structure of the topological defects. 

We first consider magnetic domain walls within one structural domain $ = 0, in which 
case Xi,2 = i>i,2- The walls separate two magnetic states related by the time reversal 
operation, so that A-ip lj2 = ±7r across the wall. Figures [1] (a-b) show two topologically 
distinct types of such walls between the degenerate B 2 states: the one in which spins in 
neighbouring layers rotate in opposite (Aipi = Aip 2 — and the same (A-0i = — Atp 2 — 
7r) directions. The free energies per unit area of the two walls are 8^S(A + C + ) and 
8-^/ S(A + C_) respectively. The domain wall with ipi = ijj 2 has a nonzero magnetic moment, 
since the magnetic configuration in the middle of the wall, (^i(O), ^2(0)) = (7r/2,7r/2), is 
of the weakly ferromagnetic A 2 type. The net magnetic moment per unit area of the wall 
is 7tMa 2 1+, where M^ 2 is the magnetisation in the A 2 phase and l± = ■sj ^ A _^ c± ^ is the 
domain-wall thickness. 

Importantly, within a domain, the reference triangles of Mn spins that are used to define 
the four magnetic phases have their apical oxygen ions tilted either towards or away from a 
common centre (Fig. [Tj). This is because the tilts of the oxygen bipyramids, described by the 
angle determine the in-plane magnetic anisotropy axes. A shift by one lattice constant 
(of the non-trimerised lattice) within the uniform domain results in a 120° rotation of the 
anisotropy axes, since $ — > $ + 2tt/3. To minimise the magnetic energy, therefore, ipi and 
?/> 2 must transform in the same way: ^1,2 4>i,2 + 2n/3. This is why the free energy Eq.([3]) 
depends on Xi,2 = i>\,2 — and why, in general, the "covariant" angles (xi, X2), rather than 
(V'ijV^)) should be used to describe the magnetic phases. 
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Now considering structural inhomogeneity, it is evident that the magnetic structure must 
respond to the presence of a structural domain wall. At a structural domain wall with 
A$ = 7r/3, for example, spins rotate either by 7r/3, in which case x is the same on both 
sides of the wall, or by — 2n/3, in which case Ax = — tt. Thus, structural domain walls 
are also magnetic domain walls. In the lowest-energy configuration, A^i = Atp2 = tt/3 
(see Fig. |4]^c)), while the next-lowest-energy configuration has Aipi = Aip2 = —2n/3 (see 
Fig. |4^d)). Domain walls in which spins in neighbouring layers rotate over different angles, 
e.g. A^i = 7r/3 and Aip 2 — —5tt/3, are higher higher in energy. It is important to stress 
that the thickness of the magnetic domain walls — l± ~ 10 2 A — greatly exceeds that of the 
structural domain wall. The A^ = — 2n/3 antiferromagnetic domain wall clamped to the 
structural boundary (Fig. |4^d)) has a nonzero moment along the c axis equal — 2tt/3Ma 2 1+ 
per unit area of the wall. 

With these considerations, we can understand the simultaneous presence of "clamped" 
and "free" antiferromagnetic domain walls in hexagonal manganites.^ Every structural do- 
main wall (A<E> = tt/3) induces a magnetic domain wall in which spins rotate by 60° or 
120° (Figs. [4] (c-d)). Because of the sign change of electric polarisation at the structural 
boundary, these antiferromagnetic domain walls appear to be clamped with the ferroelectric 
domain walls. The "free" antiferromagnetic domain walls that do not follow ferroelectric 
domain boundaries are the 180° antiferromagnetic domain walls within one structural do- 
main (Figs. [4] (a-b)). These results imply that in the minimal-energy magnetic state of the 
structural vortex spins wind around the vortex core, and the total spin rotation angle along 
a loop encircling the vortex is A0i = Aip2 = A$ = ±27r. That is, structural vortices are 
also magnetic vortices. The small-angle neutron scattering experiment on HoMnOs, the 
electric switching of magnetisation of coupled ferromagnetic/LuMnOs thin films, and the 
magnetic force microscopy study of ErMn03, all indicated the presence of an uncompensated 
ferromagnetic moment at antiferromagnetic domain walls .H2E5E2] o ur analysis summarised 
in Fig. [4] shows that both "clamped" and "free" antiferromagnetic domain walls induce mag- 
netisation along the c axis in their vicinity and several of them have a net magnetic moment. 
Although the weak ferromagnetic moment - arising from canting of Mn spins - is small, it 
can be significantly enhanced by the magnetisation of rare-earth ions coupled to Mn spins, 
as in the bulk A 2 phase 

In conclusion, multiferroic hexagonal manganites provide a rich playground for physics 
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of topological defects in multiple coexisting orders. We presented a theory based on first- 
principles calculations that explains the observed coupling between the structural distortions, 
electric polarisation and spins at the domain walls and vortices. The significance of these 
findings lies in the fact that topological defects can dominate cross-coupling responses of 
bulk materials, such as magnetoelectric switching. 

SA and MM were supported by the ZIAM Groningen under award MSC06-20 and by FOM 
grant 11PR2928. KTD acknowledges fellowship support from the International Center of 
Materials Research and computational resources provided by the CNSI Computing Facility 
at UC Santa Barbara through Hewlett-Packard. NS was supported by the ETH Zurich. 

Appendix A: First-principles calculations 

We performed ab initio calculations using Kohn-Sham density functional theory (DFT), 
as implemented in the ABINIT 45 software package J 39 | 4U | All calculations employ a supercell 
approach with periodic boundary conditions. Wave functions and charge densities are 
expanded in a plane-wave basis. Efficient computational treatment of heavy elements is 
achieved using the projector-augmented wave method for core-valence partitioning,^ which 
significantly reduces the required plane-wave energy cutoff. 

We approximate the Kohn-Sham exchange-correlation potential using the local spin den- 
sity approximation^! with a Hubbard-U correction applied to the partially-filled manganese 
d states following the Liechtenstein approach^ with double-counting corrections in the fully 
localised limit. All calculations reported here were performed with values of U — 8.0 eV 
for the Coulomb integrals and J = 0.88 eV for the intra-atomic exchange coupling, as cho- 
sen previously by Fennie and Rabe.^ We enforce an A-type antiferromagnetic ordering 
for all calculations.^ With this choice of parameters and magnetic ordering, an insulating 
electronic structure in the high-symmetry (P63/mmc) crystal structure results, with a Kohn- 
Sham band gap of 0.75 eV. Since this underestimates experimentally reported values of the 
band gapj^we carefully verified that no spurious metal-insulator transitions occur as struc- 
tural distortions are introduced, so that the free-energy landscape contains no anomalous 
features. 

The parameters given in Eq. (|2| can be extracted by considering homogeneously distorted 
periodic structures that are commensurate with the wave vectors of all distortions. The 
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smallest unit cell that can accommodate all possible values of Q, $ and V contains 30 
atoms®, or 6 formula units of YMnCV To fit the parameters, we explore the variation 
of the DFT total energy with the magnitude of Q, $ and V. The ionic positions are 
defined by the projection of the DFT ground-state P63CH1 structure onto the and K$ 
modes of P63/HIHIC combined with the desired value of Q, $ and V. However, we fully 
optimise the cell parameters for each mode distortion to eliminate stresses. The result 
is that the homogeneous contribution to stresses have been implicitly eliminated through 
renormalisation of the Landau parameters. 



The dependence of the free energy on V with Q = is shown in Fig. [2](b) , reaffirming 

the stability of the polar mode in the high-symmetry structure, and therefore the improper 
nature of the ferroelectricity in YMnC>3. Fig. [2^ a) and (c) respectively show computations 
for the trimerisation mode and the coupling between trimerisation and polar modes. 

With all homogeneous Landau parameters given in Eq. ^ specified, we now turn our 
attention to the stiffness parameters. The most convenient way to proceed is to write the 
spatial inhomogeneity in Q and V as a single harmonic, for example Q (r) = Q q e l ^ r with 
$ = 0. The stiffness energy has then the form 



This expression is harmonic in mode amplitudes Q q and V q . Note that in order to compute 
Sq, we do not need to consider spatial variations in the $ field. The computational cost of 
explicitly computing f s by imposing various short wave vectors, q, in our supercell density- 
functional calculations would be prohibitive. Instead we extract these harmonic terms us- 
ing the method of frozen phonons combined with Fourier interpolation of the inter-atomic 
force constants, a method commonly used for computing phonon band structures. Fig. [5] 
shows our calculated phonon band structure for the high-symmetry structure (P63/mmc) 
of YMnC>3. The strongest instability at the K point is the cell-tripling K% trimerisation 
mode. As previously noted, all V phonons, including the polar mode, are stable in the 
high-symmetry structure. 

Multiplying the dynamical matrix used in the computation of phonon modes by the 
weighted mass, we obtain a (f-dependent force- const ant matrix: 





(Al) 



(q) = ^M~M 3 D hJ (q) 



(A2) 
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FIG. 5: Left panel: Phonon band structure of the high-symmetry phase (P63/mmc) of YMnC>3, 
obtained using frozen-phonon methods and Fourier interpolation of the interatomic force constants. 
The most unstable zone-boundary mode at K is the K3 trimerisation mode. Right panel: Extrac- 
tion of the Sq stiffness parameter by fitting the dispersion of the force constant of the trimerisation 
mode. 

and by identifying the relevant branch, the g-dependent eigenvalues of Cy (q) are related 
to %fy or Hence, the required stiffness parameters can be extracted. As an example, 
Fig. [5] shows the extraction of Sq by fitting the q x dispersion of the force constant of the 
unstable trimerisation branch. Using this technique, we find the stiffness parameters for 
YMn0 3 listed in Table |TT} 
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